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Abstract 

A universal scaling function, describing the crossover between the Mott and 

the Efros-Shklovskii hopping regimes, is derived, using the percolation picture 

of transport in strongly localized systems. This function is agrees very well 

with experimental data. Quantitative comparison with experiment allows for 

the possible determination of the role played by polarons in the transport. 
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Electronic interactions are known to play an important role in the strongly localized 
regime. More than two decades ago, Pollak and Efros and Shklovskii pointed out 
that the long-range nature of the interactions leads to a dip in the single-particle density of 
states, p(e), at the Fermi energy. Using a constraint on the single-particle excitations, Efros 
and Shklovskii argued that this soft gap is of the form p(e) ~ e"^^^, where e is the energy 
measured from the Fermi energy and d the space dimension. The gap in the density of states 
was indeed observed by tunneling and photoemission p measurements. Moreover, it 
modifies the Mott variable-range hopping law, logi? ~ 1/T^', from x = d+ ltox = 1/2. 
The modified exponent has been observed in many experiments P]. 

The Efros-Shklovskii hopping law is expected to be relevant at low temperatures (com- 
pared to the size of the gap) and in the last few years there have been a multitude 
of experiments aimed at exploring the crossover between the Mott to the Efros-Shklovskii 
hopping regimes as a function of temperature [§]. A significant step was taken more recently, 
when Aharony and coworkers [^|10[ have argued that this crossover is described by a uni- 



versal scaling function, and obtained this scaling function phenomenologically, using energy 
additivity. It is the aim of this paper to demonstrate that the "microscopic" percolation 
picture |T1|, describing the transport in the strongly localized regime, indeed leads to such a 



universal crossover function, and to derive an equation for this function. The results agree 
excellently with existing experimental data [Q, and allow quantitative comparison between 
features in the density of the states and the transport data, taken on the same physical 
system. 

The starting point of this calculation is the mapping of the resistance problem into a per- 



colation criterion [|rT| of an equivalent random resistor network [|T^ , consisting of randomly 
placed sites. Without interactions, the activation energy from site i to site j is given by ej—et, 
where are the energy of site i, which in this case is distributed uniformly. The resistance 
between each pair of sites is given by Rij = exp {(|ej| + \ej \ + \ei — ej\)\/2T + 2rij/C,}, 
where all resistances are measured in some unit of resistance. Since the resistances vary ex- 
ponentially, the overall resistance will be determined by the weakest link, which is the largest 



resistance, i?, such that the cluster formed by all resistances (bonds), satisfying R > Rij, 
percolates. Clearly, all states participating in the percolating network (defined as occupied 
sites) must satisfy R > exp(|ej|/2T). Following |Tl]|, the percolation criterion employed here 
is the following ||T3| — given such an occupied site , the number of bonds attached to it has 
to be higher than a critical threshold, Z^, for the system to percolate, 
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2'^+ir(d/2)rf(d+l)(rf + 2)' 
= Z,T(logi?)'^+VTo, (1) 

leading directly to the Mott hopping low , logi? = (To/T)^/'^'^+^) (where the last equality 
defines Tq). In the above po is the uniform density of states, and ^ the localization length. 

In the presence of interactions, the activation energy from site i to site j is given by p| 
ej — Cj — e'^/vij, where r^j is the distance between the sites, (where all other particles are 
assumed to remain fixed). Since an activation from the ground state has to be positive, this 
leads to a constraint on the distribution of e^, and to the soft gap in the Fermi surface M. 



In fact, given that constraint, Efros [Il4| , p!5i has been able to derive a self-consistent equation 
for the density of states, 

P(ei) 1 



p(e)/po = exp dei 



(2) 



where a is some numerical constant, and an infinite band was assumed. The resulting 
density of states is clearly of the form po/(e/eo), with eo ~ {e^^poy^^'^~^\ and it can be very 
well approximated in three dimensions by f{x) = x^/(l + x^), with eo — 0.46e'^^/po. 



Knotek and Pollak [16| and Mott ITj] have pointed out that the resistance is not deter- 



mined by the single-electron density of states, but should involve many electron excitations 



(polarons). Efros [|1^] indeed demonstrated that taking into account multi-particle excita- 
tions (i.e. not assuming that all other particles remain fixed, when electron is activated 
from site i to site j), leads to an even stronger (exponential) gap in the density of states 
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1^ , [T8[] . The long-range interactions between the polarons, which are indeed the relevant ex- 
citations for the resistance, lead to a density of states which satisfies (0), and consequently 
the resistance should be of the Efros-Shklovskii type, logi? = {Tes/TY^'^. 

Taking these considerations into account, the percolation condition (Eq.(^) is now mod- 
ified in two ways. The probability of finding two sites, of energies and ej and distance 
Tij apart is proportional to p(ei)p(e2)6(e2 — ei — e^/ri2), where p(e) satisfies (H), and the 
activation energy from site i to site j is given by ej — ei — e^/vij. The resulting percolation 
condition can now be written in the form, 

I rO i'2TlogR 
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Using the one-parameter scaling form of p(e), and defining the dimensionless parameters x 
and y by 1/T = a2{e'^/^)x/el and logR = ai{e'^/C,)y/eQ, where ai are numerical amplitudes, 
this equation can be recast in the form 

' -{y'+si)/2x' Jmax{ei+{j/'-si)/2x',0} 
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with x' = a2X, y' = aiy, Si = y/y'"^ — 16a;', and S2 = y' — (£2 — Eq.(|D which is an 

implicit equation for y in terms of x, or logi? in terms of T, indeed demonstrates that the 
resistance can be written in a scaling form 

log i? = A/ (T/TJ, (5) 

in agreement with Aharony et al. 0, with A = {To/tof'^, and = {^^^ /Tof/'^ (the 
numbers were chosen such that the amplitudes of in the two asymptotic regimes 

of f{x) will agree with the amplitudes chosen by Aharony et al.). In addition to relating 
the nonuniversal amplitudes and A to eo, namely to properties of the density of states 
, which can be independently measured, Eq.(^ gives a prescription to calculating the full 
scaling function from "first principles" . 



We have calculated the crossover function in three dimensions, using Eq. and a density 
of states of the quasi-particles of the form p(e) = p^e^ + eg). In Fig. 1 we compare the 
derivative of the scaling function to the data published in Ref. [Q. In comparison we also plot 
the phenomenological function derived by Aharony et al. 0. Two points are noteworthy - the 
experimental data exhibits a rather fast crossover from the Mott regime (logi? ~ 
to the Efros-Shklovskii regime (logi? ~ The presently derived scaling function 

shows an excellent agreement with the data. The phenomenological scaling function, and 
also the scaling function derived from the procedure suggested by Lien |jlO|, show a rather 
smooth crossover, which has a poorer agreement with the data (see inset). Moreover, as one 
crosses over to the Efros-Shklovskii regime, the scaling function predicts that instead of 
a monotonic increase in the exponent from 1/4 to 1/2, there is an overshoot in the slope, a 
feature which is clearly seen in the experimental data. 

There are two physical parameters which can be extracted from the fit — Tq and eo. The 
value of Tq was found to be 1500 Kelvins, exactly as deduced in the experimental paper from 
the high temperature data. More interestingly, the best fitting value of eo was 0.04 meV, an 
order of magnitude smaller than what one would deduce from the single-electron density of 
states, measured on the same sample, by tunneling spectroscopy. This seems to be the first 
verification of the fact that it is not the single-electron density of states that determines the 
resistivity, but rather the density of states of the dressed particles, the polarons, which has 
a gap which is much smaller, due to the finite size of the polaronic cloud. The comparison 
between the present calculation and the experiment allows for the determination of the 
average cloud size, which is a few atomic lengths. 

To conclude, we have used the percolation picture to derive the crossover function from 
the Mott to the Efros-Shklovskii hopping regime. This function shows an excellent agreement 
with experimental data. Quantitative comparison with the experiment allows a quantitative 
analysis of the effect played by the many-particle excitations in the system. We hope that 
this work will stimulate more work in this direction, especially experimental investigations 
of the resistivity and the density of states in the same samples, such as those pioneered by 
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Massey and Lee 0. 

Acknowledgments: I thank Ora Entin-Wohlman for motivating me to work on this sub- 
ject, to Y. Gefen and L. Levitov for fruitful discussions and to Mark Lee for making the 
experimental data of Ref. 0] available. Work at Santa Barbara was supported in part by 
the National Science Foundation under Grant No. PHY94-07194. 



Figure Caption 

Comparison of the experimental data (diamonds) to the derivative of the scaling function 
derived in this work (continuous line), and that of Ref. |^ (broken line). The data shows 
a fast crossover, and is fitted very well by the function derived here. The inset shows 
a comparison with the previously derived functions [Ref. (broken line), and Ref. flG 
(dotted line)], both showing too smooth a crossover, compared to the data. 
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